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C^ We show that electron hopping in a lattice of molecules possessing a Berry 

phase naturally leads to pairing. Our building block is a simple molecular 

o 
o 



site model inspired by C6o ; but realized in closer similarity with Na3. In the 
resulting model electron hopping must be accompanied by orbital operators, 



Vh ' whose function is to switch on and off the Berry phase as the electron number 

changes. The effective hamiltonians (electron-rotor and electron-pseudospin) 
obtained in this way are then shown to exhibit a strong pairing phenomenon, 
by means of ID linear chain case studies. This emerges naturally from numer- 
ical studies of small JV-site rings, as well as from a BCS-like mean-field theory 
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formulation. The pairing may be explained as resulting from the exchange of 

singlet pairs of orbital excitations, and is intimately connected with the extra 

degeneracy implied by the Berry phase when the electron number is odd. The 

relevance of this model to fullerides, to other molecular superconductors, as 

well as to present and future experiments, is discussed. 
PACS numbers: 71.27.+a,74.20.Mn,74.70.Wz 
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I. INTRODUCTION 

A significant feature of the physics of unconventional superconductors such as the 
cuprates is the constraint imposed on the motion of the charge carriers by the background 
degrees of freedom, i.e. the spins in the case of the cuprate superconductors. 

In this paper we focus attention on a new class of constraints imposed on the motion of 
conduction electrons by the Berry phase ffl, or molecular Aharonov-Bohm phase ||, which 
can arise in molecular crystals with large on-site degeneracies. In general, for this to be the 
case the symmetries of electron and vibron states must match appropriately and, moreover, 
the number of electrons needs to be odd. 

As an explicit example, we have demonstrated elsewhere || the presence of a Berry 
phase in negatively charged fullerene ions Cg~ , Cg^, CgQ~ ensuing from a Dynamic Jahn- 
Teller (DJT) effect arising from coupling between the partly occupied ti u orbital and the H g 
vibron modes. 

In such a situation a physical electron (or hole) may be regarded as a composite particle, 
made up of the bare electron plus the geometrical phase which accompanies it when sitting, 
unpaired, (or more generally in a state with an odd number of electrons) on a given molecule. 
By contrast, either the absence of electrons or presence of a pair of electrons will eliminate 
the Berry phase on that molecule. We argue below that the kinematical constraints imposed 
by the Berry phase can be a factor capable of tilting the balance in favor of pairing, even 
in the presence of repulsive interactions [f§] . The way this works is for a pair of electrons on 
one molecule to gain energy by each tunnelling off onto neighbors with accompanying vibron 
excitations before coming together again. A single electron on the other hand will tunnel 
accompanied by its vibron excitations, so the electrons will gain less energy by tunnelling 
as individuals than will be gained by tunnelling as pairs. 

In this paper, we propose a class of simple coupled electron-rotor models which we believe 
capture the essential physics introduced by the Berry phase constraint. Omitting at this 
stage the complications of real Cgg" anions, our model lattice hamiltonian is instead directly 



inspired by the simpler and well known [|[ strong-coupling Berry phase molecule Na3. 

In general there are an infinite number of rotor states on each molecular site. To sim- 
plify further, we shall truncate to a 3-state model for the rotors. The resulting many-site 
Hamiltonian will be shown to take the form 

ff = -| £ AccUStSj + h.c.) (1) 

<i,j>,a 

where S + ,S~ are raising and lowering operators for a spin-1 vibron manifold for each 
molecule, and < i,j > denote first neighbors. 

To test the pairing properties of (jl]), we include a Hubbard repulsion term 

H u = U ^n^rni. (2) 

i 

We then study the half-filled state of the model numerically for small systems (4 to 8 site 
rings), and also by a mean field BCS-type calculation for the 1-D chain. 

The final conclusion is that the Berry phase coupling is found to be greatly beneficial to 
electron pairing, at least within the simple ID lattice studied so far. Pairing, in particular, 
appears to prevail and to survive even in presence of a repulsive Hubbard U, up to values 
Ufat 

In Section || we will introduce the main concepts and build our working hamiltonian, 
representing an idealized lattice of Berry phase molecules. 

Section pj is devoted to a numerical study of the 4, 6 and 8-site ID chain clusters, where 
correlations can be studied exactly, and the presence of pairing is demonstrated. 



Section |V| discusses a mean field BCS-type formulation for the infinite 1-D chain. 
We close with a discussion section, where a number of interesting open issues are also 
briefly presented. 

II. THE MODEL 

Our model system is a regular lattice of molecules. Each molecule possesses initially a 
degenerate orbital, an odd number of electrons, and a dynamical Jahn- Teller (DJT) coupling 



(with Berry phase) to a local vibration, also degenerate, of pseudorotational type. For 
simplicity, we stick to the case of one electron in a doubly degenerate orbital and one 
rotor, which provides the simplest case of Berry phase coupling. The electron can hop 
from a molecule to the next one in the usual fashion, conserving ordinary spin. We will 
generally also include an on-site ( "intra-molecular" ) electron-electron repulsion U, so that, 
for a nondegenerate level and no rotor coupling, we would have just an ordinary Hubbard 
model. Berry phase coupling to the rotors is the new ingredient giving rise to peculiar 
selection rules for electron hopping between two molecules in different rotor states. In the 
following, we describe successively the on-site hamiltonian, and the full lattice hamiltonian 
inclusive of electron hopping. 

A. On-site hamiltonian: modeling a Berry phase molecule 

We consider at each site a molecule with a partly occupied doubly degenerate electronic 
state. Suppose this orbital interacts via linear DJT coupling with a doubly degenerate 
vibration. A practical example (not relevant for superconductivity) of such a situation is 
the Na3 molecule. A single unpaired electron occupies the doubly degenerate electronic 
molecular orbital E = (E x , E y ), and the doubly degenerate vibron is a pseudorotation of the 
Na3 triangular structure. Direct spectroscopical evidence has been found ||, showing that 
the formalism we present here (supplemented by quadratic couplings which are omitted here 
for simplicity) correctly describes the dynamics of this system. 

The single-molecule, one-electron, linear coupling case, is a classic Jahn- Teller textbook 
problem ||. When only one electron occupies the degenerate state, the lowest order electron- 
vibron hamiltonian can be written: 

dq\ dql 

where the 2x2 matrices span the twofold E electronic level, and qi are the vibrational 
normal coordinates in the vibron space; the second term is the coupling between vibron and 
electron, of dimensionless strength g. 
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)+9^[ I (3) 

1/ 2 V <?2 9i 



The problem is rotationally invariant in the (91,92) space. It is therefore conveniently 
rewritten by introducing polar coordinates 9, (ft in the vibron space and a pseudospin | 
representation in the electron space: 

T\)=±^(\E x )Ti\E y )) (4) 

In terms of these quantities the hamiltonian is expressed as: 
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dq 2 qdq q 2 c 
A new total angular momentum j, defined as 



+ ^( ° *~l (S) 

9e"^ ' 



1 a /I ° \ 

j = h~ (91P2 - fcPx) - -cr, = -^ - ^ q _ J , (6) 

is found to commute with H ||. Note that a z is only a pseudospin | spanning in reality 
the twofold orbital state, and should not be confused with the true spin, which is ignored at 
this stage. Because of this pseudospin | term, the eigenvalues of j are half-odd integer, an 
amusing anomaly first pointed out by Herzberg and Longuet-Higgins [0] for the case of the 
triangular molecule. This fractionalization can be seen as a manifestation of a Berry phase 
|jlj] of 71, which the vibrons pick up from the electron degeneracy . 

Diagonalization of (|5|) must in general be done numerically ]§,§]. In the limit of strong 
coupling (g » 1), however, the massive radial 9-motion can be approximately separated 
from the 0-pseudorotation quantized by j, and both can be solved analytically in the form 
of an oscillator and a free rotor respectively. 

The resulting spectrum is classified according to j and u, the quantum number coming 
from the quantization of the radial massive motion motion: 



, s „ , L 2hu I 9 1 / „ g 8 \ 13 

E(v,j)=hu(v+-) + —lf + -±^f + ^ 1 , v = 0,1,2,. ..,j = ±-,±-,.. . (7) 

Since we shall be concerned only with the low lying rotor states, we can forget the massive 
boson ladder. Furthermore, we will express all energies in units of the pseudorotational 
quantum Q := 2*2. 



The corresponding wavefunctions in the strong coupling limit are: 

Vv,i(<7, <fi) = Vv(<?) , lw j half-odd integer (8) 

V sin 0^-2^ / 

where ip v (q) is the appropriate harmonic oscillator wavefunction, and 6j is a pseudospin 
mixing angle 
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tantf,- = 8— ^ 5* (9) 
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In the g — > oo (strong coupling) limit 0j tends to ±|. The energy can be expanded as 



9 i -.-2 , _ _i_ 4 J_ 

4 V 



E(j) = ±v+/ + 7 ±4 7I + --- ( 10 ) 



At low energies, we consider only the rotor states 



with energy: 



^ )=r *i e ^,J' j= 444- < n) 



m=j\ i = 4- 4' 4-' (12) 
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where we have omitted as usual the — ^- offset contribution (polaron energy shift), but also 
the extra zero point energy 1/4 required by the Berry phase |J. This is our simplified model 
for the one-electron Berry molecule. 

To study electron hopping among different molecules, we also need an equivalent descrip- 
tion for molecular occupancy different from one. When two electrons occupy the molecular 
orbital in a singlet state, the DJT distortion is still present. However, the orbital phases for 
the two electrons cancel each other and a Berry phase is no longer present. 

For uniformity with the one-electron case, we will still label for n=2 the pseudorotational 
levels with j. Cancellation of the geometrical phase now requires j to be integer. The 
assumed spectrum is simply that of a standard free pseudorotor: 

E(j)=f , j = 0,±l,±2,... (13) 



with wavefunctions 

1> j {<t>) = e i », j = 0,±l,±2,... (14) 

In the energy eigenvalue flT3| ) we are again omitting the polaron energy gain —\, (units of 

n). 

In the n=0 case where no electron is present, no Jahn- Teller effect and thus no pseu- 
dorotor either. However we really would like to mimic an electron-hole symmetric situation, 
like for example going Cg^ or Cqq with respect to Cg^. For this reason we assume even for 
zero electrons a pseudorotation ( |l~3l , |l4D identical to that of the two electron case. 

We finally discard occupancies higher than two; for instance, we may suppose that, if the 
average occupancy is one, states with n > 2 are strongly suppressed by Coulomb repulsions. 

To summarize, if n is the number of electrons present in the degenerate orbital, 2j 
assumes even or odd values according to whether n is even or odd. The j-dependence of the 
energy eigenvalues is quadratic, like in a free (pseudo-)rotor. The full quantum state of such 
a molecule, in the limit considered, is described by a set of three quantum numbers which 
for convenience we define as integers: 

n, (occupancy) 
m = 2j, (rotor state) 
a = 2m s , (spin state) (15) 

where m s = ±| is the z component of true electron spin. Their allowed values are constrained 
in the form: 

n = 0,l,2 

m = 21 + sin(Trn) ,1 = 0, ±1, ±2, ±3 • • • 
a = ±sin(7rn) (16) 

With reference to the physics of Cg^, discussed in Ref. 0, the similarity with the present 
model should be clear. In that case, in particular, the ground state has even L for even n, 
and odd L for odd n. 
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B. Inter-site hamiltonian — ► hopping between Berry molecules. 

To allow electrons to move among sites, we need to specify how the hopping process is 
affected by the j quantum numbers on each site. Therefore we begin considering one electron 
hopping between two neighbouring molecules. With sufficiently high point symmetry, such 
as one has for the linear chain, we will have hopping between |-E xl ) and \E Xt2 ) and between 
\Ey t i) and \E y p) only, with amplitude t x and t y respectively ||. If for instance the twofold 
degenerate state is associated to a p-orbital, then t x = t ppa , t y = t ppn , in Slater-Koster's 



notation [10]. The hopping hamiltonian then is 



Hkin — /_^ [^z( C a;,<T,l c iE,o-,2 + h.C.) + t y {c ya iCy^^ + h.C. 



(17) 



To characterize its behaviour, suppose we have on site 2 a spin up electron, with m 2 
(odd) molecular pseudospin, and we want the hopping amplitude of this electron to site 1 
with final spin up, and pseudospin mi (also odd). We need to invert relation (|J) to express 
the fermionic operators in the pseudospin basis: 



-z-v/2 

cL + cl 



y,s 



V2 



(18) 



The matrix element we compute depends also on the (even) pseudospin of the empty 
sites l 2 (final state) and l\ (initial state) in the following way: 



1 
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h 


J^kin 


T 








1 



h ni2 

t 



tx+ty / 5m 1 +l,l 1 &l2,m,2+l J r&m 1 -l,l 1 &l2,m2-'L 



+ 



tx—ty /<5m 1 +l,! 1 ^ 2 , m2 -l+<5 mi _iJ 1 (5; 2jm2 + i 

2 I 2 



(19) 



The deltas here originate through trivial orthogonality of angular momentum wavefunc- 
tions (|TT| , |T4| ), while the overall | factors originate from the a/2 factors in ( p!8|) and in (pT|). 

If we limit ourselves to the special case of intermolecular interaction with t y = t x (t y = 
—t x ), then we have an additional conservation l\ + rri2 = vn\ + l 2 (h — vni = mi — I2) of the 



total pseudospin, as shown by Eq. (|T9"D. In this cases, the spectrum separates into different 
independent manifolds. For t y = t x , they correspond to different values of total pseudospin. 
Even if somewhat unrealistic (for a real p state, t x is larger and positive, t y smaller and 
negative) this is a very convenient choice and we shall adopt it in the following. We have 
made a check, to be described in Sec. [Til], releasing this restriction, which have satisfied us 
that the physics is not fundamentally different in the more general case t y 7^ t x . Accordingly, 
we define a single effective hopping 

t = ^±^ (20) 

which in the absence of direct electron-electron interactions is the only independent param- 
eter of our model, the rotor energy quantum Q being unity. 

Similar considerations and selection rules to those discussed above for the process 
|0, 1 >— > |1, > apply to the case involving doubly occupied sites, namely |0, 2 >— > |1, 1 >, 
or |1,2 >^ |2, 1 >. 

Now we have all the ingredients to place these "Berry molecules" on a lattice to see the 
effect of local rotor coupling on electron hopping. 

C. Lattice of Berry molecules: the working hamiltonian 

As a first attempt to study this model on a lattice we consider a linear chain of N sites, 
with N e i electrons moving on them. The basis we consider is labelled by the set of all the 
rii,mi,cri, for % -- 1,---,N, so that an arbitrary state is expandable on states: 

n x n 2 ■ ■■ n N 

mi ra 2 ■ ■ ■ m^ ) ■ (21) 

0\ a 2 ■ ■ ■ a N 

These basis states are obtained by ordered applications of N e i local fermionic creation op- 
erators c] CT on a vacuum state where no electrons are present, and setting the rrii rotational 



quantum numbers to values allowed by the Berry constraint ( |16|) 
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Hopping of an electron between site i and site j is, in this space, a composite operation, 



since it implies rii —*■ m — 1, rij — > Uj + 1, but also m» — >• m» ± 1, m^ — * m.,- =F 1 ( |3~6] , p~9| 
Hence, we write a general hopping hamiltonian in second quantized language as: 



H 



kin 



"o Z^ c i,<J c j,<r{Li Lj + L i Lj 



(22) 



<i,j>,a 



where the action of the operator Lj 1 is to raise (lower) the pseudospin m^ (really an orbital 
angular momentum) by one unit: 

n x ■ 

mi 



L f 



n ; 



0"! 



rrij 



Oi 



n j 



nx ■■■ 

m>x ■ ■ ■ vfij ± 1 



(23) 



01 



O"; 



To the hopping hamiltonian H^m, we add an on-site rotor hamiltonian H rot , as well as 
an additional on-site electron-electron Hubbard interaction term Hy 

N 1 . 



Hrot "Et ^ 



3=1 



JV 

Hu = UJ2 n jA n jt[ , 
the rotor energy contribution is at site j due in the rotor state rrij being 

flrot j ~A^j ' 

as in Eq. ( jT2] , |T3"D , and having introduced a third operator 



(24) 



(25) 



nx ■ 


■ n J 


nil 


■ m 3 


ax ■ 


■ Vj 



V(Ia 



n-i 
nix 
0i 



n , 



m~ 



(J; 



(26) 



The commutation relations for these operators are [L + , L~\ = 0, [L z , L^] = ±L^. 

We stress here that the kinetic term Hj, in alone in Eq. ( [2"2"D is the relevant part of the 
new hamiltonian we want to study 

The new electron operators cj^Lf are diffrent from the original ones, cJ. jCri , c* ffi of Eq. 
(|j~7|). In particular, we have now a single band, instead of the original double band problem. 
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However, with the exclusion of all site occupancies higher than 2, all matrix elements are 
the same in the two descriptions. 

By construction, the hamiltonian (p2j ) conserves the constraints fllCf ) among the quantum 
numbers, and can therefore be diagonalized in the Hilbert space of states defined in (pT|). 
The matrix elements of the kinetic term H^m (off diagonal) and of H rot + Hjj (diagonal) on 



the basis ( pT|) are trivial, once periodic boundary conditions (PBC) are applied to indexes. 
In some cases we shall however need antiperiodic boundary conditions (ABC), replacing t 
with —t in the kinetic term involving sites 1 and N. 

The Hilbert space of the problem is infinite-dimensional even for finite N, due to the 
rrii rotor quantum numbers, which are boundless. In the numerical computations, we shall 
truncate the basis (^1|) by choosing a cutoff energy E cut , including only states having some 
local energy smaller than this E cut . Unfortunately, the choice H rot < E cut is unfair with 
respect to singly occupied states, having larger energy (at least by 1/4) than the unoccupied 
and doubly occupied lowest- j ones. To achieve better convergence even at relatively small 
E cu t, we subtract this contribution and retain those states satisfying: 

The special case E cut = is of very strong interest. Physically, this corresponds to the 
limiting case t << Q, where the intramolecular rotor energy is much larger than the hopping 
energy. In this limit, the only allowed values for rrij are (even occupancy) and ±1 (odd 
occupancy). The resulting model has six states per site, two corresponding to even (0 and 2) 
occupancies, and four to the 2x2 combinations of spin Oj and pseudospin rrij values allowed 
for one electron. It is possible and convenient to rewrite this simplified version of the model 
in terms of fictitious spin-1 states, the rrij quantum number becoming the z projection of a 
pseudospin S — 1. For this simplified version we can rewrite the hamiltonian ( |22|) replacing 
the free rotor operators L + , L~ , L z with the generators of the spin 1 algebra [TT| S + , S~ , S z 
respectively 
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H' k in=~\ E A^AStSj + h.c.) (28) 

<i,j>,a 

The full ladder of rotational states ( |TT| , |iH ) has now disappeared, being replaced just by 
the double degeneracy of the n=l sites (S z =±l), with H rot = for all n. 



The extra terms (|4]) would still need to be added to H' kin . However, H rot has the simple 
effect of giving a energy shift of Q./A per each singly occupied site. As suggested above, this 
has the same effect of an operator such as: 

ng^3.T + ^j , (29) 

i.e. it is the same as a positive Hubbard U term with U — ^, apart from a chemical 
potential. For t « Q, this term amounts to a divergent shift of the Hubbard U. Such a 
diverging term has no physical origin (the JT energy gains we have neglected in ( fLlp and 



(|T3|) are also infinite and have opposite sign). Therefore we will omit it, and simply work 
with Hamiltonian (ffij|), with the only caveat that we need to remember the shift in Hubbard 
U when comparing the results of this low-cutoff model ( p8|) with the fully converged one. 

Both hamiltonians ( Effi) and ( P8f) show a significant degree of symmetry which we can 
take advantage of. Each of them conserves the number of electrons N e i, the total pseudo- 
angular momentum 2 J = M = J2i m i (i- e - the total S z in the 5=1 pseudospin version), and 
total electron spin. For a linear chain, the lattice translational symmetry is also obvious. 
Pseudospin conservation is a result of our approximation t x = t y . The others are exact. We 
choose to study the problem ( p2[) in the manifold at half filling (N e i = N ) and at M = 
(even N), because of the higher symmetry present in this case, which includes electron- hole 
symmetry. Although we have not yet carried out a complete study of the model away from 
half filling, we believe that the basic physics will be (at least for U=0) the same, due to a 
suggestive analogy with the negative U Hubbard model which will finally emerge. 

At this point, we are set with two alternative working models. The electron-rotor (ER) 



model 



H E r = H kin + H rot + H v , (30) 
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where Hki n , H rot and Hy are given by fl22|), (|2^) , is more realistic, and is characterized by 
two parameters, the hopping energy t and the rotor energy Q. This latter quantity, in turn, 
contains the ionic mass, and will therefore make the model sensitive to isotopic changes. 
The second, electron-pseudospin one (EP) model 

H E p = Hun + H u , (31) 

where H' kin is given by fl28|) , represents the extreme molecular limit and is more idealized, 
the hopping energy t being the only parameter. Clearly, there will be no isotope effect in 
this model. 

Although the important terms are H kin and H' kin , both models are endowed with the 
Hubbard term Hu, which can describe additional repulsive interactions, and is also conve- 
nient as a gauge of the effective attractions which will arise. Having taken Q as the energy 
unit, the physical results in the ER model will depend on the two dimensionless ratios t/Q, 
U/Q. Those in the EP model will depend only upon U/t, making direct comparison with 
the simple Hubbard model particularly straightforward. 

As it will be shown, there is numerical evidence that the two models, ER and EP, lead 
to qualitatively similar effects, at least when t is not too large. Hence, it will be possible for 
many purposes to focus on the simpler EP model. 

In the next two sections, we propose to study these models on a ID linear chain, as 
follows. First, we will study numerically some very small clusters, by direct diagonalization. 
This will permit a first crude comparison between ER and EP, and also between them and 
the simple Hubbard model. Next, we will introduce a mean-field theory for model EP on the 
infinite ID linear chain. Here, the S=l pseudospin variables can be approximately integrated 
out, giving rise to negative effective electron-electron forward and backward couplings, again 
suggesting singlet pairing. 
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III. NUMERICAL STUDIES FOR SMALL LINEAR CHAIN CLUSTERS (RINGS) 

We consider here small iV-site linear chain clusters (rings), in particular iV=4,6,8, ac- 
cessible to numerical diagonalization using conventional Lanczos method. While these sizes 
are admittedly small, we find that the qualitative results for small N are clear enough at 
this initial stage [I2|. 



We define useful equal-time correlation functions for singlet superconductivity (SC) and 
spin density wave (SDW) in the standard form of g-space "structure factors": 

N 



1 

Ssc(q) = jz E ^°'-° (4l c ],T c UQ,l I»> (32) 

3,1 

1 N 

SsDw(g) = ^ E eiqU ~ l) (imii^ !-> (33) 



3,1 

We ignore alternative channels, such as charge density waves (CDW) or triplet supercon- 
ductivity ||13| , which can also be probed, but whose behaviour is not relevant at this stage. 



In particular, a CDW will be definitely disfavoured in the more general case away from half 
filling. The property of correlations ( |52| ) and (^) is that they transform into one another 
under the transformation c^f — > Cj t f, c] , —* (— l^'c^, which, remarkably, amounts simply 
to the transformation U — > —U fnj in the Hubbard model. Hence in that model Sscio) 



and SsDwiq) are perfectly symmetric around U=0, where their values must cross, the SC 
instability prevailing for U <0, the SDW for U >0. This of course is confirmed for the 
iV-site rings, as shown in Fig. |l|. In the upper panel we choose PBC, so that the free system 
(U = 0) is in a closed shell configuration of 6 electrons. In the lower panel instead ABC are 
applied. Here, the resulting shift of the single-particle fc-states yields an open shell for free 
fermions at half filling. Open shell and closed shell calculations must finally converge to the 
same answer for N —>■ oo, and their systematic comparison at finite N provides a rough but 
useful measure of finite-size corrections. 

Our strategy is therefore to calculate Ssc and Ssdw for our hamiltonian Her and Hep, 
as a function of U/t, and to find the value ofU = U* where they cross, so that superconduct- 
ing pairing prevails at all U < U*. The finding that U* is finally positive will in turn imply 
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that the bare U=0 model is approximately equivalent to a negative Hubbard U model, with 
U = U e ff, where a crude linear estimate is 

Ueff « -IT , (34) 

so long as U* is small. 

Figure |2| shows results for the EP model, obtained for iV=4,6,8 sites, at half filling 
(N e i = N) as a function of U. The two panels a) and b) correspond to the different choices 
of closed shells and open shells respectively. For instance, N = 4 and 8 correspond to closed 
shells generated with antiperiodic boundary conditions (ABC), while N = 6 does that with 
periodic boundary conditions (PBC). Conversely N = 4 and 8 yield open shells with PBC, 
N = 6 with ABC. The condensation wavevector q is correspondingly zero with PBC and 
q = § for ABC. 

These results show, strikingly, that in the EP model, a finite positive U*/t is needed to 
suppress superconductivity in favor of SDW's. Roughly, the EP model (|28| ) behaves therefore 
like a negative U Hubbard model, with U e ff oc —t For small N, the value of U e ff varies 
with N, and also depends on whether the shell is closed / open. Although we have not tried 
a systematic finite-size scaling extrapolation for U e ff to the N = oo limit, the result up to 
iV=8 suggests that 

- 0.8£ < U eff < -0.2* . (35) 

In particular at N=6 both ABC and PBC yield U e ff = — 0.37|t|, which may therefore be a 
likely value. 

We have also studied the full ER model (|3~0"D. For this model, the Hilbert space is that 



of states (|21|), with an upper cutoff in the rotor states E cut . Due to the larger Hilbert space, 
we have restricted calculations to iV = 4, 6. We proceed by calculating Ssc an d Ssdw f° r 
fixed t/Q as a function of U/fl, and we look for the value U*/Q where they cross. This again 
defines, via (|35|), a value for U e ff. Typical results are shown in Fig. [| 

Now t is an independent parameter. The effective interaction U e ff can be recalculated 
by varying t, and the results are given in Fig. [|. The main feature is that the negative 
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U e ff at small t, already found in model EP, is confirmed. Hence, model ER also leads to 
superconductive pairing for t/Q not too large. This fully confirms our expectations that 
kinematical restrictions imposed by the switching of orbital states are important in that 
regime. For larger values of t, these restrictions gradually become irrelevant, until, for 
t — > oo we recover the value U e ff — > 0. In other words, when the hopping energy is too 
large, the DJT effect does not work any more, and Fig. |4] describes how its "phase" part is 
quenched (the JT distortion magnitude is held constant in our model). 

One may suspect that the pairing we are demonstrating is just a consequence of the 
exact symmetry between the x and y degenerate molecular orbitals that we enforce by the 
assumption of having equal intermolecular hopping matrix elements t x , t y , defined in (|T7D, 
as discussed in Section |1|. A simple test releasing this assumption shows that this is actually 
not the case. If t x ^ t y , the hitherto missing term corresponding to unequal t x , t y term is 
the following kinetic additive contribution: 

H' km = ~ E 4c,, ff (L+L+ + LrL7), (36) 

<i,j>,a 

where t' is the independent hopping amplitude t' = {t x — t y )/2. This term violates the 
conservation of total m (equivalent to complete x-y symmetry). By adding to the symmetric 
hamiltonian (|30|) a term like (^) we can break continuously this symmetry, monitoring the 
effects of this on pairing, in particular on the U e ff that was defined above. 

In Fig. [5] we plot U e ff for the EP model (for simplicity) in the closed shell configurations 
for 4 and 6 sites, as a function of t/t', at fixed t + t' = t x . This figure shows that although the 
negative effective attractive U e ff is maximum for t x = t y and t x = —t y , it is not cancelled in 
the general case t x ^ t y , except for the very special case t y = (or t x = 0). The cancellation 
of pairing interaction in this limit is due to the complete breaking of the rotational symmetry, 
creating two separate bands from the x and y orbitals, the x-originated being a regular tight- 
binding band, the y one being made of localized degenerate states. Anyway, this is indeed 
a very special case, in which no t y term is present. In a more realistic situation having, say, 
t y w —t x /2, the pairing effect fully survives. 
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IV. THE INFINITE CHAIN MEAN FIELD BCS APPROACH. 

In order to get a qualitative idea of the effects of the Berry phase constraints on the 
infinite system, we first integrate out the vibron degrees of freedom in the simple spin- 
1 model Hamiltonian (|IJ). To do this, we make the further approximation of replacing 
the spin operators by pseudo-fermion operators representing the vibron excitations. We 
can then integrate out the vibron degrees of freedom and apply the BCS equations to the 
resulting interacting fermion system. 

Introducing an auxiliary spin-| fermion B ia , we classify the m — +1 vibron state as an 
"a-up" state, and the m — — 1 state as an "a-down" state. The m = state is treated as 
a vacuum state for the 5-fermions, which we will call "berryons" : 

\m = 1 >= B\\0 > , \m = -I >= B\\0 > . (37) 

Using this representation, we express the spin-1 operators by 

S+ = B \ + B i , S~ = £ T + B\ , S z = B\B^ - B\B i . (38) 

This representation is overcomplete (in particular it does not exclude unphysical states with 
Hi, rrii of different parity), and will therefore not allow a strictly variational treatment. Still, 
it is of use in exploring whether the model does or does not display tendencies toward pairing 
at the simplest mean field level. 

We rewrite the hamiltonian (|l|) in this fermion representation as: 

H ki n = ~ E 4^>(4% + 4% + B h B l + S ^T + h.c.) (39) 

<i,j>,a 

or, in Fourier representation, 

k\,k2,q,(T 

(-B 9T -Bfei-fc2+<?T + BqiB^-ki+ql + Bq^B-ki+kz-qi + ^-gl-^fei-fe+gT + h-c.) • (40) 

We take as a zeroth order mean field hamiltonian just free fermions (note that this term 
was missing in the original problem) 



h mf = Yl e kc{ a c k)CT + 2 r\ k B\ a B Ka (41) 

k,a k,a 

such that 

Emf =< Omf\H\0 M f > (42) 

is minimum. \Omf > is the direct product of a half filled Fermi sea of c-electrons and of 
a Fermi sea of berryons filled up to x = Ng/N. Here, x can be regarded as an adjustable 
variational parameter (although Hmf is not truly variational, as it violates the constraints 
required by Eq. (0)). The precise value of x is however immaterial, since the qualitative 
results we will find appear to be independent of Ng. 
Direct substitution gives 

E M f = 5-sin(-x) , (43) 

which is minimum for x = 1. The single particle excitation energies are 

e k = -26 cos k (44) 

At 

T]k = COSK , (45) 

7T 

where 6 = — sm(nx/2) is the effective mean field hopping amplitude for c-electrons. 

The next step is the determination of the first nontrivial correction to the mean field 
due to the actual interaction ( |40|) between c-electrons and berryons. These corrections are 
achieved through an expansion in the interaction (H — Hmf) around the free dynamics 
Hmf- For this purpose, we write the full many-body partition function at temperature 1/(3 



0: 



jUI ,..,.S\ \ « 7,^/,-<S>B + ki<S 2 >B-<S>] 



e- ti ) B ) c « Z MF ( e -<^>«+3l<^>B-<»>«J+-^ , (46) 



with 



ry __ rye yi* 

ZjMF — ^MF^MF 



S = f dr[H(c(r) 7 ...B(r)) - H mf (c(t), ...S(t))] 
Jo 
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<om> c = /^£*d e -/>E^4,w(flr-^)^w OMr)] 

< 0[B ka ] > B = f D ^ B ^K -So^^Ur)^^)B^(r) 0[BUT)] (4?) 

where 0[.] is any operator, and 

Z£, F = J D [clc ka ]e- to d ^ k ^UT)( dT+ e k )c kA r) (4g) 

and a similar expression for Zf /IF . 

Averaging over the non-interacting many body 5-fields in the cumulant expansion in 
( f46|) leaves an effective hamiltonian operator for the c-electrons. That expansion contains a 
first term < S > B whose form is J dr J2 c ka ( T ) c ka( r ) ■ It simply renormalizes the mean field 
parameters. The lowest order nontrivial action correction belongs to — < S 2 > B , having 
the form of an effective electron-electron interaction term 

S eff = ^ J dr J dr' Y, J2 4 lff (^) c *a<r('r) c fc 8 ff'( T ') £ **» , ( T/ )- K *i.*a.*3.*4( T ~ r ') ' ( 49 ) 

This term has a very simple significance. It corresponds to the exchange of a berryon 
particle-hole pair with singlet total pseudospin between the two electrons, as in the diagram 
of Fig. §. 

The imaginary time integration can be recast in a Matsubara frequency summation, in 

terms of a kernel 

t 2 
Kk u k 2 ,k 3 M( iuJ B) = hx+ksM+kA^Tr J2 cos ( k i + fc) cos(A;3 + k + q) 

k 
J, YJ~ 2 9k^n)9k^n + Wb) + 

g k ^u n )g ki (-uj n - u B )g ki (uj n )g k ^(-uj n + u B )} , (50) 

where q := k 2 — fci, u n are fermionic Matsubara frequencies, and g ka {u) n ) is the free fermion 
propagator in Matsubara space as defined in |15|1 . The sum over the Matsubara frequencies 



can be performed to recast Eq. (|50| ) in form 

t 2 
K klMtkzM {iu B ) = 4i+*3,ka+Jfe4 2jy : H cos (^i + k) cos(/c 3 + k + q) 

k 

[2E(w B ; £ k , £ k - q ) - E(iw B ; &, -&_,) ~ S(iw B ; -&, £ fe _ g )] , (51) 
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where T,(z,a,b) := [/f(o) — fF(b)]/[z — {b — a)}, fp() are Fermi occupation factors and 
6c : = Vk — fJ> = — ^[cos(A;) — cos(^p)] are the single-particle excitation energies for the 
berryons reduced by the corresponding chemical potential. 

We would now like to extract physical conclusions from this calculation. Since we deal 
with an effective ID electron system, we wish to use the calculated effective electron-electron 
scattering as a guide to understanding which one of the standard ID Luttinger model fixed 
points will prevail. 

In particular, for that model, an estimate of the forward and backward coupling constants 
g\ and g 2 UH w * u determine what kind of ground state to expect. 

The pair scattering amplitude we have obtained is obviously time- (or frequency-) de- 
pendent, i.e., non-hamiltonian in nature. In this sense, straight identification with true 



hamiltonian parameters such as g\ and g-i [|lt| is not automatically correct. However, we 



see no physical reason preventing us from using our derived amplitudes as effective coupling 
constants so long as we stay sufficiently close to the Fermi surface. 
We therefore identify 

g x oc K_ kF;kF!kF ,„ kF (iuj B ) (q = 7r) , g 2 oc K kFtkF _ kF _ kp (iuj B ) (q = 0) . (52) 

Direct computation of K for these special momenta gives: 

t 2 r n dk 2 
K- kF , kF , kF ,- kF (z) = -— / — sm 2 (k)[f F (-£ k ) - fp^k)}^ 

ZiV JO 7T Z^fc — Z 

K kF>kF ^ kF ,_ kF (z) = -^ [ f sin 2 (fc) [/>(-£*) - M^)]^^ (53) 

Interestingly, these two couplings arise from different terms in the hamiltonian. The for- 
mer is due to the B^B, B B* 1 terms in (|40|) , whereas the latter to the B^B^ and B B 
terms. 

In the zero frequency limit z — > 0, the two quantities are negative and coincide. In terms 
of the Luttinger model phase diagram ||13|| , this corresponds to a spin singlet superconducting 



state , which is therefore found to prevail. This is in very good agreement with the effective 
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negative U of the previous section, with the additional remark that the alternative possibility 
of charge-density waves is now explicitly ruled out. 

Actually, the zero frequency limit is singular at zero temperature, when the Fermi func- 
tions become step functions and the k integration diverges logarithmically for vanishing 
frequency around kp, the Fermi momentum for the berryons. In other words, in this ap- 
proximation gi = g2 diverge as ln|z|, at small z. This is a singular feature, due to our 
assuming the exchange of a bare, unrenormalized particle-hole pair as in Fig. |[ Higher or- 
der diagrams will modify that. More importantly, in presence of a finite pairing amplitude, 
for example, this divergence will disappear, due to a pairing gap in the berryon spectrum. 
There is in fact an exact symmetry between fermions and pseudofermions, and the two 
Cooper channels are also identical. 

V. DISCUSSION 

A very important property of our Berry phase-constrained tunnelling Hamiltonian is 
the fact that the constraint operates at the energy scale of the tunnelling matrix element 
t. Thus the pairing tendency induced by the constraints is not dictated directly by the 
strength of the intra-molecular electron-vibron coupling, but rather by the indirect effect 
of this coupling in the semi-classical limit on the relative phase space available for single 
electron tunnelling versus that for pair tunnelling. Thus the enforcing of the Berry phase 
constraint effectively separates the energy scale of the tunnelling, t, from that of the internal 
degrees of freedom of the constituent molecules. 

In the case of Cqq itself, the model is in too extreme a semi-classical limit to give a 
reasonable representation of the physics of K3C60 since, in that case, the dimensionless 
electron-vibron coupling strength g is of order ~ 0.4, whereas the strong coupling limit 
where the Berry phase representation becomes useful is for g ^ 1. 

Nevertheless, the model does illustrate a new physical principle for superconductivity in 
strongly constrained systems. It is tempting to make an analogy with the physics of the 
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t — J model of interest for describing the physics of doped Mott insulators. In that case the 
hopping Hamiltonian my be re-written as 

H t = -tJ2 {Sic\ iCjl Sj + Srct^Sf + h.c.) (54) 



where Sf, Sj' are spin raising and lowering operators. A number of recent studies ]T^| suggest 
that pairing of holes close to the half-full insulating state occurs as a result of kinematical 
constraints in this model. 

Although the present model does not have the ordered or quasi-ordered background of 
the antiferromagnetic state in the t — J model, it does have the feature that pair tunnelling 
proceeds by each partner causing a vibron excitation when executing a tunnelling step, which 
then annihilate when the pair of carriers come together again on the same site. Similarly, 
in the t — J model case, individual hole hopping is accompanied by a spin-flip, which can 
the be cancelled by the hopping of its partner. 

It is of interest to consider whether our Berry phase considerations could also apply to 
the Chevrel-phase class of superconductors such as LaMo 6 Se 8 or PbMo 6 Sj / Se 8 _, y [|T7|. In 
these materials the Mo 6 Se 8 cluster has a set of degenerate LUMO orbitals analogous to 
those in Ceo- Measurements of the doping dependence of T c indicate a sharp maximum as a 
function of doping in the unfilled LUMO shell. Thus there is the possibility of a general class 
of constraint-driven superconductors with distinctly different dependence of T c on material 
parameters than those of the conventional BCS-type electron-phonon superconductors |18] . 

Although, as stated above, we do not expect our model to be a realistic representation 
of the physics of K3C60, it would nevertheless be interesting to test experimentally whether 
the kind of electron-vibron coupling we have proposed could be observed in this compound. 
One way to do this would be through the two-vibron Raman spectrum. Our coupling 
mechanism would naturally lead to a direct electron-hole pair channel coupling to a pair of 
vibron modes. This channel would open up a gap of 2 A, in the vibron spectrum where A is 
the superconducting gap, on lowering the temperature of the material below T c . 

More generally, we observe that the pseudorotor Berry phase mechanism sketched here, 

23 



ties together electron hopping with the hopping of quanta of orbital molecular angular mo- 
mentum, which is unquenched in the free-molecule limit we start from. In the paired state 
orbital quanta are also paired, whereby the orbital excitation branch will also develop a gap 
at q = 2k,F- The gap will follow identically the superconducting gap at T = T c . The stag- 
gered orbital susceptibility should therefore be maximum at T c . In turn, the uniform q = 
orbital susceptibility may also develop a maximum, although weaker, via momentum non- 
conserving or local field effects. It is possible that orbital effects of this kind, even if weaker 
than suggested by this extreme picture, could be detectable, e.g., by NMR. In particular, 
the relaxation time 1/Ti could be enhanced at low temperature, and peak up around T c due 



to large susceptibility fluctuations |L9|. 



It would be of considerable interest to see if these new effect could be observed in K 3 C 60 , 



and RbaCgo- Encouragingly, in this latter compound, very recent NMR data |20[ seem to 
indicate a behavior of the relaxation time which is anomalous precisely in the way suggested 
above. The anomaly at T c , in particular, is seen in the Rb ion, but not on the carbon, as 
we would expect for a C 60 orbital effect. 

Finally, it is of interest to speculate that T c would be enhanced by doping our model 
system away from half-filling. Because of the nature of our pairing mechanism, the carriers 
would have more phase space for pairing if each partner in a pair could find many empty 
neighboring sites to hop on to before re-pairing. Thus the doping dependence of T c might 
be expected to have a maximum away from half-filling in systems for which this mechanism 
is driving the superconductivity. In this sense the case of half filling is probably the least 
favorable. There are indications from an exact solution for the 2-electron state that in model 
EP, U e ff is one order of magnitude more attractive near zero filling [21j . 



In the fullerides, exact half filling appears to be required by chemical stability ||22| . In the 
Chevrel systems, however, where continuous doping is feasible, one indeed finds a maximum 
of T c for a hole density close to one per molecular unit |T7fl . This corresponds to only 1/6 
filling of the narrow r 25 molecular band in that case PB[ . A second observation is that the 



correlation length should tend to be short, of the order of the intermolecular distance a, 
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since this is the scale where the energy gain takes place. In K3C60, this expectation is well 
borne out, with a correlation length of order = 2a. 

Experimentally, it would be of interest to consider building new molecular solids where 
high-symmetry Jahn- Teller molecules can exchange electrons. Larger molecules may be 
better ones because of a weaker intramolecular Coulomb repulsion U. Relatively weak JT 
coupling may provide an additional favourable circumstance, since in that case the effective 
Q is larger (although our treatment does not strictly apply there) and DJT quantum effects 
are more important. Both these conditions are met in the fullerides, but it might be possible 
to find other systems where they apply. 

VI. CONCLUSIONS 

We have proposed a model for constrained tunnelling of charge carriers in a lattice of 
Berry phase molecules, inspired by the physics of the fullerides. 

The general model (Eq. |22| ) is based on an electron-quantum rotor Hamiltonian which 
includes a (in principle infinite) manifold of vibron states on each site. Although we have 
been able to investigate the effects of a large number of vibron states on the pairing tendency, 
it is clear that the effect is strongest when only the lowest are important. For the extreme 
case where only the lowest vibron state is important (5 = 1 pseudospin model), both our 
numerical studies on small clusters (Section |I|) and our BCS-type mean field treatment 
(Section [TVD indicate a strong pairing tendency for a half-filled band. The fact that the 
model exhibits intrinsic pairing even in the presence of Hubbard repulsion U of order the 
tunnelling matrix element t is understood readily from the form of the model Hamiltonian 
(Eq. D, p2|) . On integrating out the vibron degrees of freedom, we obtain an effective BCS 
attraction of order t 2 /W, where W is a vibron bandwidth of order t. Thus the Berry phase 
constraint leads to an electron-electron attraction whose energy scale is not directly related 
to the strength of the intra-molecular electron- vibron coupling in the semi-classical limit. 

In physical terms, our model is based on the entanglement of orbital angular momentum 
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of the individual molecules with electron hopping between molecules. Pairing of electrons 
is generated by an accompanying ("singlet") pairing of orbital momenta on neighbouring 
molecules, suggesting short correlation lengths in the order of the intermolecular spacing. 
It has been argued that this mechanism might be relevant also in such other molecular 
superconductors, such as the Chevrel compounds. 
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FIGURES 
FIG. 1. Structure factors Ssc(q) ( |32[) and Ssdw(q) (p3|) for the Hubbard model (6 sites ring), 
as a function of the dimensionless parameter U/t. In the upper panel periodic boundary condi- 
tions (PBC) are applied to indexes, so that for U=0 the ground state is nondegenerate ("closed 
shell"), while in the lower panel antiperiodic boundary conditions (ABC) make the noninteracting 
ground state degenerate ("open shell"). In panel a) q = 0, while in panel b) q = tt/N, as appro- 
priate to the boundary conditions applied. Solid dots mark crossings, where the switching from 
superconductivity to spin-density waves takes place. 

FIG. 2. Structure factors Ssc (32) and Ssdw ( 133J ) for the EP model, plotted as a function of 
U/t. In panel a) we report the result for the closed shell case, while panel b) has the result for the 
open shell case, as discussed in the text. Solid dots mark crossings defining U* , indicating that 
superconductive pairing prevails even at positive U. 



FIG. 3. Structure factors Ssc ( p2| ) and Ssdw ( pq ) for the ER model (iV=4,6), plotted as a 
function of U /f2, with t = 1. In panel a) we report the result for the closed shell case, while panel 
b) has the result for the open shell case. E cu t = 6 for iV=4, and E cu t = 3 for N=6. Solid dots 
mark crossings defining U*, again indicating pairing even at positive U. 



FIG. 4. The effective Hubbard term U e ff ( pq ) for the ER model (_/V=4,6), plotted as a function 
of t/n. The iV=4 ABC and N=6 PBC are closed shells, while N=4 PBC and N=6 ABC are open 
shell cases. E cu t = 8 for iV=4, and E cu t = 4 for A^=6. For t/Vt « 1, the slope of U e ff coincides 
with U e ff/t of the EP model (Fig. |2|). For t/Q, » 1, the dynamical Jahn- Teller effects are 
suppressed, whence U e ff — ► 0. 

FIG. 5. The effective Hubbard term U e jf ( |35| ) for the EP model (A r =4,6, in closed shells) with 
the addition of the term (^), describing t x ^ t y , plotted as a function of t/t' , at fixed t+t' = t x = 1. 
The effective interaction survives everywhere except at the isolated point t/t' = 1, corresponding 
to the unphysical case t y = (or t x = 0). 
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FIG. 6. The second order effective interaction between electrons is due to the exchange of a 
pair of berryons. 
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Manini, Tosatti, and Doniach, FIG 1 
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